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ABSTRACT 

The realization that GRBs release a constant amount of energy implies that the 
post jet-break afterglow evolution would be largely universal. For a given redshift all 
afterglows should be detected up to a fixed observer angle. We estimate the observed 
magnitude and the implied detectability of orphan afterglows. We show that for rea- 
sonable limiting magnitudes (rnn m = 25) orphan afterglows would typically be detected 
from small (~ 10°) angles away from the GRB jet axis. A detected orphan afterglow 
would generally correspond to a "near-miss" of the GRB whose jet was pointing just 
slightly away from us. With our most optimistic parameters we expect that 15 orphan 
afterglows will be recorded in the SDSS and 35 transients will be recorded in a dedicated 
2m class telescope operating full time for a year in an orphan afterglow search. The rate 
is smaller by a factor of 15 for our "canonical" parameters. We show that for a given 
facility an optimal survey should be shallower, covering larger area rather than deeper. 
The limiting magnitude should not be, however, lower than ~ 23rd as in this case more 
transients from on- axis GRBs will be discovered than orphan afterglows. About 15% 
of the transients could be discovered with a second exposure of the same area provided 
that it follows after 3, 4 and 8 days for mn m = 23, 25 and 27. 



1. Introduction 

The realization that GRBs are beamed with rather narrow opening angles, while the following 
afterglow could be observed over a wider angular range, led immediately to the search for orphan 
afterglows: afterglows which are not associated with observed prompt GRB emission. Rhoads 
(1997) was the first to suggest that observations of orphan afterglows would enable us to estimate 
the opening angles and the true rate of GRBs. An expanding jet with an opening angle 9j, behaves, 
as long as its Lorentz factor 7 > Qj 1 , as if it is a part of a spherical shell (Piran, 1994). As it 
slows down and reaches 7 « 9~ the jet quickly expands laterally (Rhoads, 1999) producing a 
pronounced break in its light curve. As time progresses it can be observed over wider and wider 
observing angles, o b s - Dalai et al. (2002) have pointed out that as the post jet-break afterglow 
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light curves decay quickly, most orphan afterglows will be dim and hence undetectable. They also 
suggest that the maximal observing angle, max , of an orphan afterglow will be a constant factor 
times 6j. Hence the ratio of observed orphan afterglows, R^t^, to that of GRBs, Rqrb, will not 
tell us much about the opening angles of GRBs and the true rate of GRBs, Rqrb- 

The observation that GRBs have a constant amount of total energy (Piran et al, 2001; 
Panaitescu and Kumar 2001; Frail et al., 2001) and that the observed variability in the appar- 
ent luminosity arises mostly from variation in the jet opening angles leads to a remarkable result: 
The post jet-break afterglow light curve is universal (Granot et al., 2002). We calculate this univer- 
sal post jet-break light curve, using both first principle considerations and a calibration from the 
observed afterglows of GRB990510 (Harrison et al., 1999; Staneck et al, 1999), and GRB000926 
(Harrison et al., 2001). Using this light curve we estimate the maximal flux at an observing angle 
9 bs from the jet axis. Using this flux we estimate the total number of orphan afterglows that can be 
observed given a limiting magnitude and the distribution of these orphan afterglows as a function 
of o b s and the redshift, z. 

The assumption that the total energy is constant implies that orphan afterglows will be de- 
tected roughly up to a constant observing angle 9 max (which is independent of 9j, for 9j < 9 max ). 
In this case R°^ ph oc Rqrb anc ^ can teach us about the distribution of 9j (Granot et al. 2002). 

We describe our analytic model and the phenomenological fits to the observations in §2. We 
present our estimates for the observed rate of orphan afterglows in §3. We stress that these are 
idealized optimistic estimates that do not consider various observational obstacles. We do not 
consider how and whether these transients could be identified as afterglows and distinguished from 
other transients. In §4 we compare our estimates for the expected rate of orphan afterglows with 
the capabilities of several surveys. 



2. The Model 

Consider an adiabatic jet with a total energy E and an initial opening angle 9j. We consider 
a simple hydrodynamic model for the jet evolution (Rhoads 1999; Sari, Piran & Halpern, 1999). 
Initially the jet propagates as if it were spherical with an equivalent isotropic energy of 2E/9 2 : 

E = (2ir/3)9 2 R 3 i 2 nm p c 2 , (1) 

where n is the ambient number density. The spherical phase continues as long as 7 > 9J 1 . At this 
stage the jet expands sideways relativistically in the local frame, with 7 = _1 and adiabaticity 
implies that 

E = (2Tr/3)R 3 nm p c 2 , (2) 
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and the radius of the shock remains constant . Note that the evolution during this phase is 
independent of 6j, whose only role is to determine the time of the jet break. If E and n do not 
vary from one burst to another then the light curve during this phase will be universal, depending 
only on the microscopic parameters (the equipartition parameters, €B, e , and the power law index 
p of the electrons distribution) of the specific afterglow. During both phases the observed time is 
given by: 

t = (1 + z)R/Acj 2 . (3) 
Equations 2 and 3 yield that the jet break transition takes place at (Sari, Piran & Halpern, 1999): 

t jet = 0.7(1 + z)(£ 5 i/no) 1/3 (0,/O.l) 2 days . (4) 

where Q x denotes the value of the quantity Q in units of W x times its (c.g.s) units. Due to relativistic 
beaming, an observer located at o b s outside the initial opening angle of the jet (9 bs > 0j) will 
(practically) observe the afterglow emission only at tg when 7 = 0~ h l g : 

t e = A(e obs /6j) 2 t jet . (5) 

Roughly at tg the emission is also maximal. From then on it decays in the same way as for an 
on-axis observer. The factor A in Eq. 5 is of order unity. The value of A is uncertain and it differs 
from one model to another. For example, according to model 2 of Granot et. al (2002) A < 1 (i.e. 
the peak flux seen by observer at o i, s , is seen before the jet opening angle is 8 b s ), while A > 1 
according to their model 1. Therefore, lacking a better knowledge we use A = 1 throughout this 
paper. 

The synchrotron slow cooling light curve for the initial (spherical) phase is given by Sari, Piran 
& Narayan, (1998) and modified by Granot & Sari (2002). Sari Piran & Halpern (1999) provide 
temporal scalings for the maximal flux and the synchrotron and cooling frequencies during the 
modified hydrodynamic evolution after the jet break. Combining both results [using the Granot k, 
Sari (2002) normalization for the fluxes and typical frequencies] we obtain the universal post jet- 
break light curve. For the optical band which is usually above the typical synchrotron frequency, 
but can be either above or below the cooling frequency, we find: 

F v>Vc {t) = 460 [so(p)/so(2.2)] lO^e^eJ^ 

n^E^Z 7 t2 ys ^l{ 2 (l + zf^Dll % M Jy , (6) 
F u<Uc (t) = 170[ 5l (p)M(2.2)] 10— <l\e B *_ 2 

n^E^ 7 t2 v A7 )/2 {l + z)^Dll ^Jy , (7) 



1 More detailed calculations show that R increases slowly and 7 decreases exponentially with R (Rhoads 1999; 
Piran 2000) 
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where Dl is the luminosity distance and go(p) = lCT a56p (j> — 0.98) [(p — 2)/(p — l)] p_1 , gi(p) = 
10~°- 31p (p — 0.04) [{p — 2)/{p — The cooling frequency, for which equation 6 and 7 give the 

same flux is: 

p_o. 98 y p)/2 

p - 0.04 J 

e-fA^aftl + ^^Hz. (8) 

Note that z^ c does not depends on 0j, and it remains constant in time (after the break). 
Using Eq. 5 we obtain the maximal flux that an observer at 6 b s will detect: 

F^ZiOots) = 1670 [ Sb (p)/ flb (2.2)] A-v?~\if_ 2 

3p— 2 

n^El&v u f(l + z?-*l*DZ*fi^_ x ^Jy , (9) 

^Pota) = 620 [ ffl (p)M(2.2)] A-^\ef_ 2 

nfE 50 . 7 v[lj»\l + z)^D^fi^_ x ^Jy • (10) 

One notices here a very strong dependence on 6 \, s . The peak flux drops quickly when the observer 
moves away from the axis. Note also that this maximal flux is independent of the opening angle of 
the jet, 9j. 

Once we know F™ ax (0 o b s ) we can estimate, 6 max (z,m), the maximal observing angle at which 
an afterglow is brighter than a limiting magnitude, m. We then proceed and estimate the total 
number of orphan afterglows brighter than m. 

We can estimate the observed flux using the "cannonical" values of the parameters. However, 
these are uncertain. Alternatively, we can use the observations of the afterglows of GRB990510 
(Harrison et al., 1999; Stanek et al., 1999) and of GRB000926 (Harrision et al., 2001), both showing 
clear jet breaks to obtain an "observational" calibration of F u (0 o b s ). As the 9 b s dependence in 
Equations 9 and 10 is similar we can write: 

F(9 obs ) = F f(z)9-f , (11) 

where -Fo is a constant and f(z) = (1 + z) l+ P Djj^% includes all the cosmological effects. j3 is the 
spectral slope {F v oc v@). The value of j3 depends on whether the observed frequency is below or 
above the cooling frequency v c : (3 u<Uc = (1 — p)/2 and (3 u>Uc = —p/2. For simplicity, we consider 
throughout the paper only the case where v > v c and (3 = —p/2. 

If all GRBs are similar with the same total energy, ambient density and microscopic parameters 
(e e> B andp) than Fo is a "universal" constant. However, its value is uncertain. We need the observed 
flux and at a certain time to determine Fo. We use the on-axis observed afterglows at the break 
time, when an observer at o b s = 6jA~ 1 / 2 (for a narrower jet) would observe a similar flux to an 
on-axis observer (for the observed jet). Using the parameters for GRB990510 and GRB000926 we 



v c = 3.62 x 10 15 
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estimate^ directly. For GRB990510 (Harrison et. al. 1999) we have 9j = 0.08, p = 2.1, z = 1.6 
and rribreak = 20mag, where mb rea k is the R band magnitude of the afterglow at the break time 
and F (990510) = 0.01 fiJy. For GRB000926 (Harrison et. al. 2001) 9j = 0.1, p = 2.2, z = 2 
and mbreak = 20mag and Fo(000926) = 0.03 fiJy. Both values are rather close to the theoretical 
estimate with E 5 i = 0.5, n = 1, e e = 0.1, e# = 0.01, p = 2.2 which yields F = 0.02 [iJy and 
magnitude 20 at the break time (for a burst at z = 2 with 9j = 0.1). 

Granot et al. (2002) show, using a more refined simulation of the off-axis light curve (their 
model II), that when the off-axis light curve 'joins' the on axis light curve, the off-axis maximal flux 
is a factor of few less than the on-axis flux. For example a burst at z = 1 with £51 = 0.5, n = 1, 
e e = 0.1, 6b = 0.01, p = 2.5 is estimated in this model at 24th magnitude at 9 b. s = 0.22. This 
result corresponds to Fq = 0.002 fi Jy. According to the observations and the fact that the off-axis 
maximal flux is a factor of few less than the on-axis flux we use in the following a "canonical" model 
of Fo = 0.003 jJL Jy and p = 2.2. We should keep in mind that there is a large uncertainty (factor 
of ~ 10) in the absolute value of the flux. This uncertainty does not change 9 max that scales as 
F^ 2p by a large factor. However, the overall detection rate depends, almost linearly on Fq. The 
total number of observed bursts depends not only on ^ aa . oc F^ p but also on the duration, t b s 
(see Eq. 12 below) which is oc 6^ aa ,, Together we find that the rate is oc F^ p . 

For a given limiting magnitude, m, we can calculate now the total number of observed orphan 
afterglows. For a given z we define 9 max (z,m) such that 23.6 — 2.5Logio{F u [9 max (z,m)]} = m 
where F u is the observed flux in units of fiJy. Figure 1 depicts contour lines of the inverse function 
m(6 b s , z) on the (6, z) plane. We use throughout this paper a "standard" cosmological model with 
n M = 0.3, ^ A = 0.7 and h = 0.7. 

According to our model an observer at 9 f, s will observe the orphan afterglow if 9 Q b s < 
8max(z,m). The afterglow will be detected over t b s (z,9,m): 

At ■ * 

t obs (z,9,m)*-^(9l ax -9* obs ) (12) 

until the signal decays below the limiting magnitude. In this model which we denote model A, 
an orphan afterglow could be observed from a solid angle 2 2ir[9 max (z, m) 2 — 9?] where 9 max (z,m) 
does not depends on 9j. This model, and more specifically eq. 11, are quite generally valid for 
9 max > 9j. However if 9 D b s ~ 9j, FJJ iax {9 b s ) is more sensitive to the exact jet model being used. 
Specifically, in such a case there will be a significant contributions to the signal at 9 a b s even before 
the jet's opening angle reaches 9 Q b s . This leads to a larger flux at 9 Q b s and hence to a larger 9 max . 
For example, Granot et al. (2002) find, for their model II, that the solid angle between 9j and the 
real maximal angle (9 max ) is approximately constant even when 9 Q b s exceeds the asymptotic value 
of 9 max (given by equation 11), i.e. 9^^ = 9^% + 9?. To take care of this effect we consider also 
model B, in which we assume that this relation holds. While model A is more conservative, model 



2 The factor of 2 reflects the expectation that the jets are double sided. 
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Fig. 1. — Left panel: R magnitude contour lines (from 17- left bottom to 29) for Fq = 0.003 /^Jy. 
The right panel depicts the GRB rate n(z)(dV/dz)/(l + z) as a function of z (for z pea k = 1), as an 
indication to where we should expect most GRBs. As we plot only the region with z > 0.5 the flux 
depends almost as a power law on the redshift. 
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B allows for more orphan afterglows. We see later that for mn m > 25 surveys, both models give 
comparable results. Model B gives a factor of 2.5 more orphan afterglows for mn m = 23 (with our 
"canonical" parameters). It also predicts a larger ratio of orphan to on axis afterglows. 



where n(z) is the rate of GRBs per unit volume and unit proper time and dV(z) is the differential 
volume element at redshift z. We assume that the GRB rate is proportional to the SFR and is given 
by: n(z) = B 10°- 75z for z < z pea k and n(z) = B 10°- 75z P eafe for z pea ^ < z < 10. The normalization 
factor, B, is found by the condition: Rq^ b = fbRcRB = Jo° {dV / dz)n(z) / {I + z)dz where /& is 
the beaming factor and Rqrb = 667 yr -1 . In the following we consider two star formation rates 
peaking at z pea k = 1 and at z pea k = 2 (Bouwenn, 2002). 

Usually the detector's exposure time is smaller than t o b s (z,0,m). Thus, the number of de- 
tectable orphan afterglows in a single snapshot, over the whole sky, is: 



3. 



Results 



The rate of observed orphan afterglows (over the entire sky) is (in model A) : 




(13) 




(14) 



If we consider now model B we find: 




(15) 



where Fu m is the limiting flux for detection. In this case R or phl 'R-grb ~ (Ro/Riim) 1 ^ p is independent 
of fb. Using t bs(9 bs) ~ te(0 max ) — tg(0 o b s ), we estimate the number of observed orphan afterglows: 




(16) 



We note that since F oc A-p , both and N (») scale as 1/A 
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The z integrand of Eq. 14 gives the normalized redshift distribution of the observed orphan 
afterglows for model A (see Figure 2a). This distribution has a broad flat region all the way from 
z = to z pea k . It sharply declines above z pea k ■ For a SFR with a higher z pea k there are significantly 
fewer orphan afterglows (see Figure 3). The function peaks (weakly) at z pea k with higher limiting 
magnitude as more orphan afterglows are observed in this case. For model B the results are similar. 

The 6 integrand of Eq. 14 [with 9 max (z 1 m) replaced by z max (6,m) and integration performed 
over z\ yields the angular distribution for model A (see Figure 2b). For 25th limiting magnitude 
the median observing angle is 5°-7° depending on the SFR rate. With lower z pea k most GRBs are 
nearer and hence stronger and detectable to a slightly larger angles. These values are close to the 
jet opening angles. Most of the observed orphan afterglows with this limiting magnitude will be 
"near misses" of on-axis events. The O max (z pea k) are significantly larger (10° and 14° respectively). 
With 27th magnitude the median of the angular distribution moves way out to 12°. This is larger 
than most GRB beaming angles but still narrow, corresponding only to 2% of the sky. Again the 
results for model B are similar (as long as 6j < 0.1). 

As the sky coverage of GRB detectors is not continuously Air it is possible that a GRB pointing 
towards us (i.e. with o b s < 9j) will be missed, but its on-axis afterglow will be detected. In Table I 
we provide the ratio of the number of orphan afterglows to total number of afterglows (both on-axis 
and orphan), in a single snapshot of the sky. As expected this ratio is large for small jet opening 
angles and for large limiting magnitudes (and models A and B give similar results in this limit) and 
it decreases with increasing 9j and decreasing mu m . An insensitive search (mn m = 23) yields for 
model A, a comparable (or even larger) probability that a visual transient is a missed on-axis GRB 
whose afterglow is detected or an orphan one. However, for model B, this ratio is less sensitive to 
either 0j or mi im for 9j ~ max , and even for mi im = 23 the majority of the afterglows would be 
orphans. 

Figure 3 depicts the number of orphan afterglows per square degree (and in the entire sky), in 
a single exposure, as a function of the limiting magnitude. The thick lines are for model A (with 
various parameters) while the thin lines are for model B. In the "optimistic" case (Fq = 0.015 /iJy, 
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Table 1: Table I: The ratio of the number of orphan afterglows to the total number of afterglows 
(both on-axis and orphan), in a single snapshot, for different limiting magnitudes and jet opening 
angles. The first value is for model A, while the second one in brackets is for model B. 
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9j = 0.05) the predictions of both models are similar. With the "canonical" normalization (Fq = 
0.003 /iJy, 2 pea fc = 1 and 9j = 0.1) model B predicts 2.5 times more orphan afterglows than model 
A, for a limiting magnitude of 23. These numbers provide an upper limit to the rate in which 
orphan afterglows will be recorded as pointdike optical transients in any exposure with a given 
limiting magnitude. These numbers should be considered as upper limits, as our estimates do not 
include effects such as dust extinction that could make these transients undetectable. 

We ask now what will be the optimal strategy for a given observational facility: short and 
shallow exposures that cover a larger solid angle or long and deep ones over a smaller area. The 
exposure time that is required in order to reach a given limiting flux, Fi im , is proportional to F^. 
We can now divide the number density of observed orphan afterglows (shown in Figure 3) by this 
time factor and obtain the rate per square degree per hour of observational facility. Figure 4 depicts 
this rate with a calibration to a hypothetical 2m class telescope that reaches a 25th magnitude with 
one hour exposure. Other telescopes will have another normalization but the trend seen in this 
figure remains. Figure 4 shows that shallow surveys that cover a large area are preferred over deep 
ones (in both models). This result can be understood as follows. Multiplying eq. 16 by Ff im shows 
that the rate per square degree per hour of observational facility oc Ff in ^^ p . As p > 1 (2 — 2/p > 0) 
a shallow survey is preferred. A practical advantage of this strategy is that it would be easier 
to carry out follow up observations with a large telescope on a brighter transient detected in a 
shallow survey. Additionally one can expect that the number of spurious transients will be smaller 
in a shallower survey. However, if one wishes to detect orphan afterglows, one should still keep the 
surveys at a limiting magnitude of ~ 23, as for smaller magnitudes the number of on-axis transients 
detected will be comparable to (model B) or even larger than (model A) the number of orphan 
afterglows (see Table I). 

The number density of orphan afterglow transients, given by Eq. 14 and depicted in Figure 3 
gives an indication of the number of single events recorded above a limiting magnitude. Clearly, 
a single detection is not enough in order to identify that the transient indeed corresponded to an 
afterglow. A second following detection of the transient with a decrease by, say, 1 magnitude, would 
increase the probability that the transient is an orphan afterglow. We estimate the rate of a double 
detection of this nature by an afterglow search with a limiting magnitude mu m and a time delay 
dTobs between the two exposures of the same region in the sky. We consider, therefore, a detection 
of a transient in one exposure and a second detection after time dT Q b s in which the transient has 
decayed by more than one magnitude but is still above the limiting magnitude of the survey. This is 
of course an idealized situation and we do not consider realistic observational problems (like weather 
etc.) which may make the detection rate lower. Not surprisingly, the ratio of a double detection to 
a single detection depends on the time delay, dT b. s (see Figure 5). The optimal time delay is the 
time in which an afterglow from z pea k and with mu m decays by one magnitude: dT b s = 3 days for 
mn m = 23, dT b s = 4 days for mn rn = 25 and dT^g = 8 days for mu m = 27. With this optimal 
choice, the fraction of the detected transients that will be detected in a second exposure, is 20% for 
mu m = 27 and 12% for mu m = 23. Note however, that the curve becomes narrower for low limiting 
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magnitudes, making the exact timing of the second exposure more critical. Deeper surveys are less 
sensitive to the choice of the time delay. 

It is interesting to compare the rate of observed orphan afterglows to the true GRB rate as 
a function of the beaming factor, /& « 2/0? . At this point models A and B differ the most. For 
model B this ratio is oc (Fo/Fn m ) 1 / p and for the "canonical" parameters it holds the values of 
0.004, 0.014, 0.031 for mi im = 23, 25, 27 independently of /;,. In model A this ratio vanishes for low 
fb values. For high /& values the ratio converges to model's B constant ratio. For the "canonical" 
parameters this ratio in model A reaches half of the asymptotic value for ff, > 50, 125, 250 for 
mu m = 27, 25, 23 respectively. Thus, for strong beaming (or high mu m ) it is possible to estimate 
the true GRB rate from a determination of the orphan afterglow rate. For a low beaming factor 
this ratio is model dependent. 

4. Discussion 

We have calculated the number per unit solid angle of orphan afterglows that could be detected 
by idealized surveys with different limiting magnitudes. Our calculations are based on a simple 
idealized model for the hydrodynamics of the sideways expanding jets. Light curves from other 
models, including numerical simulations (Granot et al., 2001) described in Granot et al., (2002) show 
similar behavior. As the models differ when o b s ~ 6j we also consider an alternative assumption 
for this range. This yields analytic expressions for R orp h and N orp h (equations 15 and 16). We have 
shown that both models yield similar results for surveys deeper than mu m = 25 and have minor 
differences for surveys with mu m ~ 23. 

The normalization factor, Fq, of this light curve (Equation 11) is somewhat sensitive to the 
choice of the model parameters. Even observed afterglows with a clear jet break do not yield a well 
defined normalization factor, Fo = A~ p f(z)~ 1 9^ p F L ,(tj et ,9 i )S = 0), because of the uncertainties in 
9j, A, and F v (tj et ,9 ohs = 0). As N orph oc AF^ P oc F v {tj et ,9 ohs = 0) 2 / p /A (equation 16), we expect 
the uncertainty in N orp h to be similar to that in Fo (a factor ~ 10). However, this uncertainty does 
not influence the trends that indicate a strategy for optimal orphan afterglow survey. 

We stress that we do not consider here practical observational issues, such as dust extinction 
or the ability to identify a weak transient on top of a background host galaxy. We have also 
assumed idealized weather conditions: when considering a "verified" identification by two successive 
detections we assumed that a second exposure is always possible. All these factors will reduce 
the observed rate by some unknown factors below our optimistic values. We do not consider how 
recorded transients could be identified as afterglows. For example when we discuss a single detection 
we mean a single photometric record of a transient in a single snapshot. This is of course a very 
weak requirement and much more (at least a second detection) would be needed to identify this 
event as an afterglow. 
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Finally one should be aware of the possible confusion between orphan afterglows and other 
transients. Most transients could be completely unrelated (such as supernovae or AGNs) and 
could be distinguished from afterglows with followup observations. However, there could be a class 
of related transients like on-axis jets whose gamma-ray emission wasn't observed due to lack of 
coverage or failed GRBs (Huang et al. 2002). The only quantitative measure that we provide 
is for the ratio of on-axis afterglows to truly orphan ones (i.e. off-axis afterglows whose prompt 
gamma-ray emission would have been detectable if viewed from within the jet opening angle). 

We have shown that for a given dedicated facility the optimal survey strategy will be to perform 
many shallow snapshots. However, those should not be too shallow as we will observe more on-axis 
afterglows than off-axis ones. A reasonable limiting magnitudes will be ~ 23 — 25. With this 
magnitude we should perform a repeated snapshot of the same region after 3 to 4 days. 

We consider now several possible surveys. The SDSS (York, 2000) has mn m ~ 23 over 10 4 
square degrees in the north galactic cap. Figure 3 suggests that under the most optimistic condition 
we expect 15 detection of orphan afterglow transients in the north galactic cap of the SDSS. A 
careful comparison of the SDSS sky coverage and the exposure of relevant GRB satellites could 
reveal the rate of coincident detection (GRB - Sloan optical transient) from which one could get a 
better handle (using Table I) on the possible rate of detection of orphan afterglows by SDSS. Under 
the more realistic assumptions we expect a single orphan afterglow transient in the SDSS. Note 
that the 5 filters SDSS system could possibly select orphan afterglow candidates even with only 
a single detection (Rhoads 2001; Vanden Berk et.al, 2001), though this method could be applied 
only to bright sources. As the north galactic cap survey provides only a single snapshot with no 
spectral information at magnitude of 23, even if a transient will be recorded in this survey it will 
be impossible to identify it as a transient. Clearly this is not the best survey for orphan afterglows 
search. However, The south galactic cap survey, observes the same part of the sky repeatedly with 
mn m ~ 23. It can achieve a higher magnitude, 25th, for steady sources. If the repetition interval 
is longer than 2 weeks, than for most orphan afterglows, the second observations can be considered 
as a new one. Therefore, given that the time spent in the south cap is comparable to the time 
spent in the north cup, the detection rate should be comparable. The Advantage of the south cap 
is the ability to identify transients sources even if the decay process could not be observed, and the 
possibility of ruling out sources that are not transient in nature, such as AGNs. 

Consider now a dedicated 2m class telescope with an aperture 0.5D° of mi im = 24 for a 10 
minutes exposure and mu rn = 25 for a 1 hour exposure. Under our most optimistic assumptions 
it will record up to 35 orphan afterglows per year in the shallow mode and a third of that (13 
afterglows) in the deeper mode. Using our "canonical" model we find two orphan afterglows per 
year in the shallow mode. 

The VIMOS camera at the VLT has a comparable aperture of 0.06D but it can reach mn m = 26 
in 10 minutes and mu m = 27 in an hour. A single mi im = 27 frame would have a 0.002 orphan 
afterglows under the most optimistic assumptions (1.1 • 10~ 4 canonically) . Thus, at best 2-3 orphan 
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afterglows should be hiding in every 1000 exposures. The Subaru/Suprime-cum has a similar 
sensitivity but its field of view is four times larger. Therefore every 1000 exposures should contain 
no more than 10 orphan afterglows.lt will be however a remarkable challenge to pick up these 
transients and confirm their nature from all other data gathered in these 1000 exposures. 

We thank Tom Broadhurst, Re'em Sari and Avishay Gal- Yam for helpful remarks. This 
research was supported by a US-Israel BSF grant and by NSF grant PHY-0070928. 
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Fig. 2. — Top panel (a): the normalized redshift distribution of observed orphan afterglows in 
a single snapshot (thick lines) and the integrated z distribution (thin lines). Bottom panel (b): 
the angular distribution of observed orphan afterglows in a single snapshot (thick lines) and the 
integrated o t, s distribution (thin lines). The circles depict 9 max (z pea k 1 m) and the squares depict 
(0 o bs). In both panels we use Fq = 0.003 /uJy, and 6j = 0.1. The different curves correspond to 
miim = 25, z peak = 1 (solid line); mu m = 25, z peak = 2 (dashed line); m Um = 27, z peak = 1 
(dashed- dotted line) and mn m = 27, z peak = 2 (dotted line). 
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Fig. 3. — The number of observed orphan afterglows per square degree (left vertical scale) and in 
the entire sky (right vertical scale), in a single exposure, as a function of the limiting magnitude 
for detection. The thick lines are for model A with three different sets of parameters: i) Our 
"canonical" normalization Fq = 0.003 ^Jy, z pea k = 1, 0j = 0.1 (solid line). The gray area around 
this line corresponds an uncertainty by a factor of 5 in this normalization, ii) Our most optimistic 
model with a relatively small 0j = 0.05 and a large Fq = 0.015 /xJy (dashed- dotted line), hi) The 
same as our "canonical" model, except for z pea k = 2 (dotted line). The thin lines are for model B, 
where the solid line is for our "optimistic" parameters, while the dashed line is for our "canonical" 
parameters. Both models are similar for the "optimistic" parameters while model B predicts slightly 
more orphan afterglows then model A for the "canonical" parameters. 
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Fig. 4. — The detection rate of orphan afterglows per square degree per hour of telescope time. 
The curves that correspond to the "canonical" normalization are for a single detection in model A 
(solid line), a single detection in model B (dashed line) and a double detection in model A (dotted 
line). The rate is calibrated for a telescope that reaches mu m = 25 in a one hour exposure. It can 
be scaled up or down for other telescopes. 
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Fig. 5. — The ratio of double detection of orphan afterglows to a single detection, as a function of 
the time delay, dT b s , between the two exposures. The three models are: Fq = 0.003 /iJy, z pea k = 1 
and mu m = 23 (dashed line) mn m = 25 (solid line) and mn m = 27 (dashed-doted line). 



